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Abstract 

The ground state structures of the A^^i^xC wurtzite (WZ) ahoys with x =0.25, 0.5, and 0.75 
are revealed by a ground state search using the valence-force held model and density-functional 
theory total energy calculations. It is shown that the ground state WZ alloy always has a lower 
strain energy and formation enthalpy than the corresponding zinc-blende (ZB) alloy. Therefore, 
we propose that the WZ phase can be stabilized through alloying. This novel idea is supported by 
the fact that the WZ AlPo.sSbo.s, AlPo.75Sbo.251 ZnSo.5Teo.5, and ZnSo.75Teo.25 alloys in the lowest 
energy structures are more stable than the corresponding ZB alloys. To our best knowledge, this 
is the first example where the alloy adopts a structure distinct from both parent phases. 

PACS numbers: 61.50.Ah,61.66.Dk,64.70.kg,71.15.Nc 
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III-V and II- VI semiconductors usually crystallize into one of two forms: hexagonal 
wurtzite (WZ) and cubic zinc blende (ZB) structures. The ZB and WZ structures have the 
same local tetrahedral environment and start to differ only in their third-nearest-neighbor 
atomic arrangement. Despite the structural similarity, there are some significant differences 
in the electronic and optical properties Compared to the hexagonal structure, the 

cubic phase has a more isotropic property, higher carrier mobility, lower phonon scattering, 
and often better doping efficiency. In contrast, the WZ phase has a larger band gap (usually 
direct), a spontaneous electric polarization, and a lower propagating speed of dislocations 
and thus an improved lifetime of the laser diodes 3|]. For certain device applications, one 
phase is preferred over the other. To have a controllable way to synthesize the desired phase, 
it is important to understand the mechanism for stabilizing a certain structure. 

In general, the WZ structure is preferred over the ZB structure when the ionicity of 
a compound is high ^. This is because the ideal WZ structure has a larger Coulomb 
interaction energy with a larger Madelung constant, whereas the ZB structure leads to 
a better covalent bond formation sl, Q, 8]. To change the stability, one often grows 



materials into different forms. For example, many ZB compounds can adopt the hexagonal 



WZ structure when forming nanowires (NWs) 
that the stability of the WZ NW is due to the 



St 



111 ]. Empirical calculations suggested 



than the ZB NW with a similar diameter 



12 



act that the WZ NW has less surface atoms 



13j. Theoretical calculations also showed 



that stability of WZ compounds such as GaN can be changed when carriers are introduced 



through dopin; 



ing 



Moreover, metastable phases can be synthesized by employing non- 



equilibrium growth techniques. For example, metastable ZB GaN can be grown on cubic 



substrates 



If 



In this paper, we show for the first time that the ground state (GS) WZ alloy (WZA) 
always has a lower strain energy than the corresponding ZB alloy (ZBA). Therefore, if strain 
energy is dominant in alloy formation, stable GS ternary WZAs can form even though the 
binary constituents are more stable in the ZB phase. This provides an opportunity to form 
desired WZAs through alloying. Our first principles calculations confirm this idea, showing 
that WZ AlPo.sSbo.s, AlPo.75Sbo.25, ZnSo.5Teo.5, and ZnSo.75Teo.25 have lower total energies 
than the ZB counterparts. 



The GS structures of ZBAs have been extensively studied 
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2lj. For 



instance, it was shown that the GS ZB A0.5B0.5C alloy (Without loss of generality, B ion 
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is assumed to have a larger radius than A ion, and C could be anion or cation) adopts the 



tetragonal chalcopyrite structure (space group I42d, No. 122) [l7|. However, the knowledge 



22 1 showed that the 



of the GS structures of WZAs remains incomplete. Our previous work 
GS structure of the A0.5B0.5C WZA is of the /9— NaFe02 type with the space group Pna2i 
(No. 33) as shown in Fig. [U^a). Here, in this work, we identify that the GS structures of 
A0.25B0.75C and A0.75B0.25C WZAs have the structures shown in Fig. [T](c) or (d) with the 
space group P2i (No. 4). 

The formation enthalpy of isovalent semiconductor alloys Aa.Bi_a;C is defined as 

AHf = E{x) - [xEac + (1 - x)EbcI (1) 

where Eac-, Ebc, and E{x) are the total energies of bulk AC and BC, and the A^jBi-^-C alloy 
with the same crystal structure (WZ or ZB). It is well known that for lattice- mismatched 
isovalent semiconductor alloys, the major contribution to the formation enthalpy is the 



strain energy. The strain energy (Eg) could be described well by the VFF model (231. l24l. |25| . 
which considers the deviation of the nearest-neighbor bond lengths and bond angles from 
the ideal bulk values. Here, we consider all possible supercells with up to 32 atoms per 
unit cell. For each supercell, we consider all possible configurations of alloys with x=0.25 
and 0.75. The VFF model is used to relax the structure and predict the energy of the 
configuration. We considered Ga^Ini-^N, AlP^Sbi^^., ZnS^Tei-a;, and GaPx-^si-x- They 
have various degrees of lattice mismatch: 10.1%, 11.5%, 11.8%, and 3.7%. Our calculations 
reveal the Lazarevicite structure (space group Pmn2i, No. 31) shown in Fig. [U^b) has the 
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lowest strain energy for A0.25B0.75C for all four different sets of VFF parameters 
The Pmn2i A0.25B0.75C structure has the same supercell as the Pna2i A0.5B0.5C structure. 
One can get the Pmn2i A0.25B0.75C structure by replacing one half of the A atoms in the 
Pna2i A0.5B0.5C structure with B atoms so that each C atom has one neighbor A atom and 
three neighbor B atoms. The Pmn2i WZ A0.25B0.75C structure is similar to the famatinite 



ZB A0.25B0.75C structure 



20] in that they both have similar local environment for C atoms. 
For the WZ A0.75B0.25C alloy, we identify two low strain energy structures with the P2i 
space group (No. 4) [P2i-I: Fig. [D^c), and P2i-II: Fig. [11(d)] . In contrast to the Pmn2i 
A0.75B0.25C structure, there are some C atoms which have four A neighbor atoms in both 
structures. In this sense, the F2i-I and P2i-II WZ A0.75B0.25C structures are similar to the 



HQ 



Q8 and Q16 ZB A0.75B0.25C structures [17|, [l9(]. As shown in Tabled the P2i-I structure 



has the lowest strain energy. However, the strain energy difference between the P2i-II and 
P2i-I structures is very small, less than 0.3 meV/atom. 

To see if the GS structures predicted by VFF strain energy calculations are consistent to 



the density 
tions 
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unctional theory (DFT) total energy calculations, we performed DFT calcula- 
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30| on the WZ A0.25B0.75C and A0.75B0.25C alloys with the Pmn2i, P2i-I, 



and P2i-II structures. Our results are shown in Tabled We can see that the Pmn2i struc- 
ture is not the GS of the WZ A0.25B0.75C alloy because the P2i structures have a slightly 
lower total energy, even though the Pmn2i structure has a lower strain energy. This can 
be explained in terms of the Coulomb interaction. For the A^.Bi^i.C alloy, the charge of 
A ions is different from that of B ions due to the different electronegativity. In this case, 
the Coulomb interaction is found to stablize the P2i structures over the Pmn2i structure 

n 

because the P2i structures has larger charge fluctuation [31]. Similar situation also occurs 
in ZBAs [21]. 

After knowing the GS structures, we now compare the strain energy of the ZBA and 
WZA using the VFF model. Our results are shown in Table [Tll We can see that for all 
considered systems (Ga^^Ini-^^N, AlPa.Sbi_2^, ZnSa;Tei_a;, and GaPi^Asi-^. with x = 0.25, 
0.5, and 0.75), the GS WZAs always have a lower strain energy than the GS ZBAs. The 
difference in the formation enthalpy mainly depends on the size of the lattice mismatch of 
alloy: For the first three Ax.Bi_^.C alloys with large lattice mismatch (Aa > 10%), the strain 
energy difference dEg at x = 0.5 is around 5 meV/atom, whereas, the difference dEg for 
GaPo.sAso.s (Aa < 4%), is only 0.7 meV/atom. 

Our above VFF calculations show that the WZ structure has a better ability to accomo- 
date the strain in a lattice mismatched alloy than the ZB structure. This is due to the fact 
that the WZ structure has a larger degree of freedom to release the strain. First, for the 
binary compound, the four-atoms unit-cell WZ structure has three free parameters (a, c, u). 
In contrast, the two-atoms unit-cell ZB structure only has one free parameter (a). Second, 
the WZA is also more flexible than the ZBA. As an example, we compare the 16-atoms 
WZ Pna2i and 8-atoms ZB chalcopyrite structures. In both structures, each C atom bonds 
with two A and two B atoms. In the ZB A0.5B0.5C chalcopyrite structure, there are three 
free parameters. However, there are fifteen free parameters in the WZ Pna2i structure. 
The larger number of degree of freedom in the WZ Pna2i structure leads to an enhanced 
flexibility in strain relaxation. 
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For a better understanding of the strain relaxation in WZAs, we can also decompose 
the total strain energy into the contributions from each atom [32!]. In this way, we can 
tell which kind of atoms are mainly responsible for the different behavior between the WZ 
and ZB alloys. This analysis shows that the main difference comes from the B ions with 
a large size. For example, the total contributions to the strain energy in the chalcopyrite 
(Pna2i) AlPo.sSbo.s alloy (here A=P, B=Sb, and C=A1) from Al, P, and Sb are 23.7 (22.6) 
meV/atom, 2.2 (1.4) meV/atom, and 7.4 (1.9) meV/atom, respectively. We can see that the 
strain energy difference from Sb ions contributes 74% to the total strain energy difference. 
In addition, we find that the difference mainly comes from the deviation of the Al-Sb-Al 
bond angles from the ideal value (109.47°). In chalcopyrite AlPp.sSbo.s alloy, the maximum 
deviation of the Al-Sb-Al bond angles is 5.4°, much larger than that (2.7°) in WZ AlPg.sSbo.s 
alloy. 

The calculated DFT formation enthalpy difference dAHf = AHf{WZA) - AHf{ZBA), 
where AHf{WZA) [AHf{ZBA)] is the formation enthalpy of the WZA (ZBA) defined in 
Eq. [H are shown in Table [Tll We see that it follows the same trend as the strain energy 
difference, i.e., the GS WZA always has lower formation enthalpy than the corresponding 
ZBA. However, the lower formation enthalpy in the WZA does not necessarily mean that 
the WZA has lower total energy than the ZBA because the formation enthalpy are defined 
with respect to the pure bulk compounds with the same lattice structure, whereas the total 
energy difference between the WZ and ZB A^Bi^^C alloys should also include the bond 
energy difference (dEb) between the WZ and ZB phases of the parent binary compounds. 
We define Ewz-zb{AC) [Ewz-zb{BC)] as the energy difference between the WZ and ZB 
phases of the AC (BC) compound. The bond energy difference dEi,{x) between the WZA 
and ZBA as a function of x are then defined as: 

dEhix) = xEwz-zb{AC) + (1 - x)Ewz-zb{BC). (2) 

The total energy difference between the WZA and ZBA can then be calculated as 

dEtot = dAHf + dEb (3) 

It is clear from Eq. (3) that only when the formation enthalpy difference (dAHf) is more 
negative than —dEb, the WZA can be more stable than the ZBA. 

The DFT total energy calculations are performed to determine which alloy structure is 
the GS phase of Gaa;Ini_a;N, AlPa;Sbi_a;, ZnSa;Tei_a;, and GaPa;Asi_3; with x = 0.25, 0.5, 
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and 0.75. For the parent compounds, we find that the energy differences Ewz-zb between 
the WZ and ZB phases are -5.6, -10.8, 3.5, 6.5, 3.2, 6.0, 8.8, and 11.4 meV/atom for GaN, 
InN, AlP, AlSb, ZnS, ZnTe, GaP, and GaAs, respectively. In agreement with previous first 
principles calculations [5] and experimental observations, we find that GaN and InN have 
the WZ GS structure, whereas the other compounds take the ZB phase as the most stable 
structure. The DFT results from the alloy calculations are summarized in Table [Tll For 
alloys with WZ binary constituents (InN and GaN) or small lattice-mismatched ZB binary 
constituents (GaP and GaAs), the GS alloy structure (Ga^Ini_2-N and GaP^jAsi-^^) is the 
same as the parent compounds. However, for AlPo.sSbo.s, AlPo.75Sbo.25, ZnSo.5Teo.5, and 
ZnSo.75Teo.25, the WZA structure is the GS phase despite that the alloys are formed from 
ZB parent compounds. It is interesting to note that compounds such as MnTe (CdO), which 



has the stable NiAs 
ZB compounds [33 



Rocksalt) structure can be stabilized in the ZB phase by alloying it with 



34( 1 . Here we show that the alloy can be stabilized in a structure that 



is different from both parent structures. This remarkable alloy stabilized wurtzite structures 
originate from the fact that the gain in the strain energy relaxation when forming the WZA 
is larger than the average of the bond energy difference between the ZB and WZ phases. 
For example, dAHf = —6.50 meV/atom and dE^ = 5.01 meV/atom for AlPo.5Sbo.5. It is 
also interesting to see that, the alloy stabilization energy dAHf for A0.75B0.25C is larger than 
A0.25B0.75C, i.e., the WZA is more favored when a large atom is mixed into a smaller host 
than a smaller atom is mixed into a large host. 

In order to determine the concentration x at which dEtot < 0, the dependence of the 
difference in the formation enthalpy [dAHf{x)] between the WZA and ZBA on the concen- 
tration X is essential. By definition, dAHf{0) = and dAHf{l) = 0. The x dependence 
of dAHf{x) can be obtained by fitting the data in Table. [TTlto a fourth order polynomial. 
The fitted result for the AlPj.Sbi_a; alloy is shown in Fig. [2l We can see that the curve is 
asymmetric with respect to x = 0.5; the minimum of dAHf{x) occurs at x = 0.61. Following 
Eq. [31 we obtain the dependence of dEtot on x (Fig. [2]). It is seen that the minimum of dEtot 
occurs at X = 0.66. And when 0.34 < x < 0.89, the WZ AlPi.Sbi_K alloy is more stable than 
the ZBA. For ZnSi.Tei_a., the result is similar, and the lowest concentration and highest 
concentration for a stable WZ ZnS^Tei_a; alloy are 0.39 and 0.87, respectively. 

In summary, we have identified the GS structures of the A^Bi_^,C WZAs with x = 0.25, 
0.5, and 0.75. Using VFF and DFT calculations, we show that the GS WZA always has 
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a lower strain energy and formation enthalpy than the corresponding ZBA, and thus the 
strain relaxation favors the formation of the WZA. We confirm this idea by showing that 
GS WZ AlP^Sbi_^ (ZnS^Tei_^) with 0.34 <x< 0.89 (0.39 <x< 0.87) is more stable than 
the corresponding ZBA although their parent structures crystallize in the ZB phase. 
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DE-AC36-08GO28308. The work in Fudan (FU) is partially supported by the National 
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Special Funds for Major State Basic Research, and Postgraduate Innovation Fund of FU. 



7 



[1] C.-Y. Yeh, S.-H. Wei, and A. Zunger, Phys. Rev. B 50, 2715 (1994). 

[2] M. van Schlfgaarde, A. Sher, and A.-B. Chen, J. Cryst. Growth 178, 8 (1997). 

[3] L. Sugiura, J. Appl. Phys. 81, 1633 (1997). 

[4] A. Garcia and M. L. Cohen, Phys. Rev. B 47, 4215 (1993). 

[5] Chin-Yu Yeh, Z. W. Lu, S. Proyen, and A. Zunger, Phys. Rev. B 45, 12130 (1992); ibid 46, 
10086 (1992). 

[6] J. C. Philhps, Bonds and Bands in Semiconductors (Academic, New York, 1973). 

[7] J. St. John and A. N. Bloch, Phys. Rev. Lett. 33, 1095 (1974) 

[8] J. R. Chehkowsky and J. C. Philhps, Phys. Rev. B 17, 2453 (1978). 

[9] G. Patriarche, F. Glas, M. Tchernycheva, C. Sartel, L. Largeau, and J.-C. Harmand, Nano 
Lett. 8, 1638 (2008). 

[10] M. Koguchi, H. Kakibayashi, M. Yazawa, K. Hiruma, and T. Katsuyama, Jpn. J. Appl. Phys. 
31, 2061 (1992). 

[11] C. X. Shan, Z. Liu, X. T. Zhang, C. C. Wong, and S. K. Hark, Nanotechnology 17, 5561 
(2006). 

[12] T. Akiyama, K. Sano, K. Nakamura, and T. Ito, Jpn. J. Appl. Phys. 45, L275 (2006). 

[13] V. G. Dubrovskii and N. V. Sibirev, Phys. Rev. B 77, 035414 (2008). 

[14] G. M. Dalpian and S.-H. Wei, Phys. Rev. Lett. 93, 216401 (2004). 

[15] G. M. Dalpian, Y. Yan, and S.-H. Wei, Appl. Phys. Lett. 89, 011907 (2006). 

[16] V. K. Lazarov, J. Zimmerman, S. H. Cheung, L. Li, M. Weinert, and M. Gajdardziska- 

Josifovska, Phys. Rev. Lett. 94, 216101 (2005). 
[17] S.-H. Wei, L. G. Fcrreira, and A. Zunger, Phys. Rev. B 41, 8240 (1990). 
[18] L. G. Ferreira, S.-H. Wei, and A. Zunger, Phys. Rev. B 40, 3197 (1989). 
[19] Z.W. Lu, D. B. Laks, S.-H.Wei, and A. Zunger, Phys. Rev. B 50, 6642 (1994). 
[20] J. Z. Liu, G. Trimarchi, and A. Zunger, Phys. Rev. Lett. 99, 145501 (2007). 
[21] S. Chen, X. G. Gong, and S.-H. Wei, Phys. Rev. B 77, 073305 (2008). 
[22] H. J. Xiang, S.-H. Wei, J. L. F. Da Silva, and J. Li, Phys. Rev. B 78, 193301 (2008). 
[23] P. Keating, Phys. Rev. 145, 637 (1966). 
[24] R. Martin, Phys. Rev. B 1, 4005 (1970). 



8 



[25] J. L. Martins and A. Zunger, Phys. Rev. B 30, 6217 (1984). 

[26] K. Kim, W. R. L. Lambrecht, and B. Segall, Phys. Rev. B 53, 16310 (1996). 

[27] Our first-principles density functional theory (DFT) calculations were performed on the basis 



of the projector augmented wave method 



2811 encoded in the Vienna ab initio simulation 



package [29| using the local density approximation For relaxed structures, the atomic 

forces are less than 0.01 eV/A. 
[28] P. E. Blochl, Phys. Rev. B 50, 17953 (1994); G. Kresse and D. Joubert, ibid 59, 1758 (1999). 
[29] G. Kresse and J. Purthmiiller, Comput. Mater. Sci. 6, 15 (1996); Phys. Rev. B 54, 11169 

(1996). 

[30] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

[31] R. Magri, S.-H. Wei, and A. Zunger, Phys. Rev. B 42, 11388 (1990). 

[32] For each bond involving the atom i, half of bond stretching energy is contributed to the atom 
i. The bond bending energy contributions from atom i include all bond angles which are 
centered at atom i. 

[33] S.-H. Wei and A. Zunger, Phys. Rev. Lett. 56, 2391 (1986). 

[34] Y. Z. Zhu, G. D. Chen, H. Ye, A. Walsh, C. Y. Moon, and S.-H. Wei, Phys. Rev. B 77, 245209 
(2008). 



9 



FIG. 1: (a) The GS Pna2i structure of the WZ A0.5B0.5C alloy, (b) The Pmn2i structure, which 
is the lowest strain energy structure of the WZ A0.25B0.75C alloy, (c) The lowest strain energy 
structure (-P2i-I) of the WZ A0.75B0.25C alloy, (d) The low strain energy structure (P2i-II) of the 
WZ A0.75B0.25C alloy. 
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FIG. 2: (Color online) Differences in the formation enthalpy (dAHf), total energy (dEtot), and 
bond energy (dEt,) between the WZ and ZB AlPajSbi-a; alloys. 
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TABLE I: VFF-calculated strain energy (in meV/atom) of WZ Ga^Ini_j;N, AlP^^Sbi-a;, 
ZnS^iTei-^;, and GaPa;Asi_a; alloys for the Pmn2i, P2i-I, and P2i-II structures at x = 0.25 and 
0.75. The numbers in parenthesis are the DFT calculated formation enthalpies. * and | indicate 
the GS structures obtained from the VFF and DFT calculations, respectively. 

Structures Pmn2i P2i-I P2i-II 

Gao.25Ino.75N 12.54* (13.13) 14.13 (12.38^) 14.17 (12.50) 
Gao.75Ino.25N 20.16 (16.91) 19.01* (13.25*) 19.28 (13.47) 
AlPo.25Sbo.75 19.43* (18.88) 21.02 (18.19*) 21.20 (18.37) 
AlPo.75Sbo.25 24.41 (27.09) 23.46* (23.35) 23.69 (23.32*) 
ZnSo.25Teo.75 13.19* (21.36) 14.19 (20.76*) 14.35 (21.05) 
ZnSo.75Teo.25 15.14 (29.36) 14.73* (26.96) 14.81 (26.86*) 
GaPo.25Aso.75 2.32* (2.26) 2.42 (1.99) 2.44 (1.96*) 
GaPo.75Aso.25 2.57 (2.48) 2.55* (2.09*) 2.58 (2.12) 
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TABLE II: Differences in the VFF strain energy (dEs), DFT formation enthalpy (dAHf), DFT 
bond energy (dEh), and DFT total energy {dEtot) between the GS WZAs and ZBAs. Energy is in 
meV/atom. 

dEs dAHf dEb dEtot 
Gao.25Ino.75N -4.03 -5.41 -9.53 -14.94 

Gao.5Ino.5N -4.61 -5.79 -8.21 -14.00 
Gao.75Ino.25N -4.77 -7.25 -6.91 -14.16 
AlPo.25Sbo.75 -6.18 -4.68 5.76 1.08 

AlPo.5Sbo.5 -7.35 -6.50 5.01 -1.49 
AlPo.75Sbo.25 -5.52 -6.14 4.27 -1.87 
ZnSo.25Teo.75 -4.42 -3.89 5.33 1.44 

ZnSo.5Teo.5 -5.18 -5.43 4.62 -0.81 
ZnSo.75To().2o -3.G7 -5.12 3.90 -1.22 
GaPo.25Aso.75 -0.62 -0.78 10.75 9.97 

GaPo.5Aso.5 -0.74 -1.05 10.12 9.07 
GaPo.75Aso.25 -0.39 -0.85 9.48 8.63 
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